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The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope provides both direct and indirect 
measurements of galactic cosmic rays (CR). The LAT high-statistics observations of the 7 GeV - 1 TeV 
electron plus positron spectrum and limits on spatial anisotropy constrain models for this cosmic-ray 
component. On a galactic scale, the LAT observations indicate that cosmic-ray sources may be more plen- 
tiful in the outer Galaxy than expected or that the scale height of the cosmic-ray diffusive halo is larger 
than conventional models. Production of cosmic rays in supernova remnants (SNR) is supported by the 
LAT y-ray studies of several of these, both young SNR and those interacting with molecular clouds. 

© 201 1 Elsevier B.V. All rights reserved. 


1. Introduction 

Among the earliest stimuli for studying cosmic y rays was the 
recognition by Hayakawa [1 ] and Morrison [2] that cosmic ray par- 
ticles interacting with interstellar matter would inevitably produce 
high-energy y rays. The first high-energy y-ray sky map by OSO-3 
[3] confirmed this expectation. Gamma-ray and cosmic-ray stud- 
ies, both observational and theoretical, have been closely linked 
ever since, with principal results for photon energies above 
30 MeV coming from satellite detectors (SAS-2 [4], COS-B [5], 
EGRET/Compton Gamma Ray Observatory [6]) and ground-based 
TeV telescopes (for a recent review, see [7]). 

While sharing many experimental techniques, cosmic-ray and 
y-ray studies offer complementary approaches. Cosmic-ray detec- 
tors measure the particles directly in the solar neighborhood, 
including detailed spectral and composition analysis. Observations 
with y-ray telescopes study cosmic-ray interactions at a distance, 
providing indirect insight about the particles. The y-ray observa- 
tions have the potential to trace sources of cosmic rays, since the 
photons are not deflected by galactic magnetic fields as the 
charged particles are. 

The present review concentrates on results from the Large Area 
Telescope on the Fermi Gamma-ray Space Telescope satellite, which 
was launched on 11 June 2008. The Fermi LAT offers three lines of 
research related to cosmic rays: 
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1 . Because the LAT is a particle detector, it can make direct mea- 
surements of cosmic rays, especially cosmic-ray electrons 
(CRE). 

2. The large-scale y-ray properties of our own and other galaxies 
offer some insight into the distribution of cosmic rays and their 
sources. 

3. The LAT measures the y-ray properties of likely acceleration/ 
interaction sites of cosmic rays such as supernova remnants. 

2. The Fermi large area telescope 

The LAT is a y-ray telescope [8], designed to distinguish photons 
in the energy range 20 MeV to more than 300 GeV from the high 
background of energetic charged particles at the 565 km altitude 
orbit of the Fermi satellite. The tracking section has 36 layers of sil- 
icon strip detectors to record the tracks of charged particles, inter- 
leaved with 16 layers of tungsten foil (12 thin layers, 0.03 radiation 
length, at the top or Front of the instrument, followed by 4 thick 
layers, 0.18 radiation length, in the Back section) to promote y- 
ray conversion into electron-positron pairs. Beneath the tracker 
is a calorimeter comprised of an 8-layer array of Csl crystals 
(1.08 radiation length per layer) to determine the y-ray energy. 
The tracker is surrounded by segmented charged-particle anticoin- 
cidence detectors (ACD; plastic scintillators with photomultiplier 
tubes) to reject cosmic-ray background events. The LAT’s high sen- 
sitivity stems from a large peak effective area (~8000 cm 2 ), wide 
field of view (~2.4 sr), good background rejection, superior angular 
resolution (68% containment angle ~0.6° at 1 GeV for the Front sec- 
tion and about a factor of 2 larger for the Back section), and excel- 
lent observing efficiency (keeping the sky in the field of view with 
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scanning observations). In scanning mode, the LAT observes the 
entire sky every 2 orbits (about 3 h). Pre-launch predictions of 
the instrument performance are described in [8], 

In addition to being a y-ray telescope, the LAT is intrinsically an 
electron detector. One of the biggest analysis challenges for elec- 
tron studies is the separation of CRE from the much more abundant 
hadrons (primarily protons). Key to the background rejection is the 
LAT capability to discriminate between electromagnetic and ha- 
dronic showers, based on the different event topology in the three 
subsystems. This is largely in common with the photon analysis, 
the most prominent difference being the use of the ACD. 

Unlike previous generations of pair conversion telescopes, the 
LAT basically triggers on all the charged particles crossing the active 
volume, with no harsh photon selection built into the hardware trig- 
ger. This drastic change of approach is largely due to the use of silicon 
detectors, allowing precise tracking with essentially no detector-in- 
duced dead time. Some event filtering is necessary prior to downlink 
in order to fit the data volume into the allocated telemetry band- 
width. The onboard filter is really tailored to preserve the y-ray effi- 
ciency by discarding events that are unlikely to be y rays. The rate of 
events depositing more than 20 GeV into the calorimeter, however, 
is small enough that we can afford effectively to disengage the filter 
above that energy, transmitting to the ground high-energy charged 
particle events along with the y rays. 

3. Direct cosmic ray measurements 

A century after the pioneering cosmic-ray studies by Domenico 
Pacini and Victor Francis Hess we have a huge body of knowledge 
about these energetic particles. In the LAT energy range, many 
instruments have measured the energy spectra of CR particles 
and their composition, including not only protons and electrons, 
but also heavier nuclei. They can be roughly divided into inclusive 
measurements (i.e. p and e~ + e + spectra) and measurements of the 
charge composition (i.e. p/(p + p) and e + /(e~ + e + )), providing com- 
plementary information on the origin, acceleration and propaga- 
tion of cosmic rays. This dichotomy has a parallel in the variety 
of experimental techniques routinely exploited since the mid- 
1960’s, in that cosmic-ray instruments generally fall in one of the 
two categories: calorimetric experiments (that cannot distinguish, 
per se, the charge sign) and magnetic spectrometers. 

Fig. 1 is a partial compilation of such instruments, classified by 
acceptance 1 and effective observation time. It is not a trivial task to 
give an exhaustive and unbiased account of the huge experimental 
effort that turned the first observations into precision science. In 
many cases the geometry factor is energy dependent, and some 
modern instruments (particularly the magnetic spectrometers) can 
operate in different configurations with different acceptances. In 
addition, each detector is sensitive in a different energy band (for 
any given observable), which makes a fair comparison extremely dif- 
ficult. Despite these caveats, Fig. 1 is accurate enough to illustrate 
some interesting aspects of the historical development of cosmic- 
ray measurements: calorimetric experiments feature, in general, a 
larger acceptance than magnetic spectrometers; and space experi- 
ments have the obvious advantage of a relatively longer collection 
time. Fig. 1 also puts Fermi into the context of this development: 
the LAT’s unique combination of large acceptance and long observa- 
tion time results in an unrivaled exposure factor 2 exceeding 10 s 
m 2 sr s, as indicated by the diagonal gray lines. 


1 The acceptance, or geometry factor (C F ), is the integral of the effective area over 
the instrument field of view and constitutes the proportionality factor between the 
event rate in the detector and the incident flux. 

2 The exposure factor is the product of the geometry factor and the effective 
observation time and is directly related to the total number of events recorded by a 
given instrument. 


4. Direct CR measurements with Fermi 

High-energy cosmic-ray electrons and positrons constitute a 
peculiar component of the cosmic radiation in that, unlike protons 
and heavier nuclei, they rapidly lose energy by synchrotron radia- 
tion on galactic magnetic fields and by inverse Compton (IC) scat- 
tering on the interstellar radiation field. In the energy range of 
interest for the LAT— i.e. above a few GeV, where losses dominate 
over escape— the average lifetime for electrons against synchrotron 
and IC losses is r e (E) oc £~\ whereas the corresponding escape time 
for protons scales with the inverse of the diffusion coefficient: 
t p (£) oc E s (with 5 a; 0.3-0.6). This implies that the all-electron 
spectrum at Earth is significantly steeper than the proton spec- 
trum. More importantly, at energies above a few hundred GeV 
the majority of the electron and positron flux must be contributed 
by sources closer than a few hundred pc. In other words, high-en- 
ergy CRE really probe cosmic-ray production and propagation in 
nearby or local galactic space. 

The other piece of conventional wisdom is that unlike electrons, 
which are accelerated in CR sources, the vast majority of positrons 
arise as secondary products of cosmic-ray interactions in our Gal- 
axy. In the absence of primary sources of positrons it is straightfor- 
ward to demonstrate that in the simplest scenarios the positron 
fraction is expected to decrease with energy like £~ <s , e.g. [9]. As 
we shall see in a moment, this expectation is at odds with the most 
recent measurements [10]. 

Fermi began its operations at about the same time two impor- 
tant measurements [11,10] were published that challenge the sim- 
plest models of cosmic-ray production and propagation. These 
reports stirred an enormous excitement in the community, mainly 
in connection with the possibility of detecting signs of new physics 
in the context of indirect dark matter searches. In 2009 the Fermi 
LAT collaboration published the first systematics-limited measure- 
ment of the inclusive spectrum of cosmic-ray electrons and posi- 
trons [12] between 20 GeV and 1 TeV. Though harder than 
previously thought, the Fermi spectrum does not show evidence 
of any prominent feature, such as the one reported by the ATIC 
experiment [11]. This conclusion was later independently con- 
firmed by the H.E.S.S. [13] and Pamela [14] experiments, as well 
as by a dedicated analysis of the Fermi data aimed at the highest 
end of the energy range [15], Fig. 2 summarizes the current status 
of measurements. 

The Fermi measurement— especially when the low-energy exten- 
sion of the spectrum published in [15] is taken into account— is in 
tension with a single power-law model, though this cannot be defi- 
nitely ruled out within the systematic uncertainties. This tension 
has motivated a number of papers in which a consistent interpreta- 
tion of the Fermi inclusive spectrum and the rise in the positron frac- 
tion above r; 10 GeV measured by Pamela [10] is attempted through 
the introduction of an additional nearby source of high-energy elec- 
trons and positrons (for a more detailed discussion see [17] and ref- 
erences therein). It is fair to say that, to date, the issue is not settled, 
and there is no firm consensus on whether observations really indi- 
cate the existence or nature of this extra component. Indeed there 
are alternative scenarios (see for instance [18,19]) that might ex- 
plain the observations without invoking any new (astrophysical or 
exotic) source of high-energy electrons and positrons. 

For the future, with more than 5000 electrons above 1 TeV per 
year crossing the LAT, Fermi can hope to provide an independent 
confirmation of the spectral cutoff measured by H.E.S.S. [20] in this 
energy range. Even with a selection efficiency on the order of 10% 
(in line with the figures of merit of the current analysis), the statis- 
tics will not be an issue. At the same time, the collaboration is ac- 
tively working on the energy measurement, which is the most 
challenging part of the analysis. 
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Fig. 1. Compilation of the characteristics (in terms of geometric factor and collection time) of some balloon and satellite cosmic-ray and y-ray detectors. Care has been taken 
in extracting the relevant information from the original publications in a consistent fashion. We believe that the information in the plot is accurate within a factor of «2, 
which is adequate for the illustrative purpose of the plot itself. Filled markers indicate space-based experiments, while open markers indicate balloon experiments; 
conversely squares indicate calorimetric experiments, while circles indicate magnetic spectrometers. 



Fig. 2. Cosmic-ray electron plus positron spectrum measured by Fermi in 1 year of observation, adapted from [15]. Some of the most recent measurements are shown for 
comparison. The dashed line is a pr e-Fermi model, shown for reference. Many of the data points are taken from the cosmic-ray database maintained by Strong and 
Moskalenko [16]. 


Without an onboard magnet, Fermi LAT cannot directly distin- 
guish the charge sign. However the Earth’s magnetic field creates 
natural shadow regions from which particular charges are forbid- 
den— simply because their paths are blocked by the Earth— and 
can be effectively exploited in order to separate positrons from 
electrons. The idea of using the geomagnetic field to study the 
charge asymmetry in the leptonic component of cosmic rays is 
not new. In fact the possibility of using the east-west asymmetry 
in the geomagnetic cutoff rigidity has been explored in several bal- 
loon experiments (see [21] and references therein for more infor- 
mation). In addition to its large field of view and its excellent 
directional capabilities, Fermi can take advantage of the accurate 
models of the Earth’s magnetic field [22] and the particle tracing 
computer programs that are nowadays readily available. One of 
the main challenges in the analysis is the fact that, as the energy 
increases, the shadow regions become smaller and smaller. Cou- 
pled with the falling spectrum, this effect in practice limits the 


highest accessible energy to «200 GeV. In addition to that, these 
regions are located near the Earth’s horizon, which is outside or 
at the limits of the LAT field of view in the nominal sky-survey 
operations. Preliminary results of this analysis are in agreement 
with the positron fraction measured by Pamela and will be the sub- 
ject of a forthcoming paper. With the AMS detector now installed 
on the International Space Station, new data in the unexplored 
high-energy range will be available in the near future. 

The electron component of cosmic rays is not fully described by 
the inclusive energy spectrum and the positron fraction: the de- 
gree of anisotropy in the arrival directions is an additional key ob- 
servable, especially in a scenario in which the inclusive electron 
spectrum is dominated by one (or a few) nearby sources. Some de- 
gree of anisotropy is indeed expected even in a plain diffusive sce- 
nario in which the sources are continuously distributed, just due to 
the fact that the Earth is not at the center of the diffusion halo. It is 
natural to expect that the flux anisotropy can be much larger if a 
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Fig. 3. Fermi LAT 95% CL upper limits for the integral dipole anisotropy S as a function of the minimum energy (adapted from [23]). For comparison the expected anisotropy 
levels for two pulsars (Vela and Monogem) that are plausible candidate sources of high-energy electrons are shown. The values of the model parameters are described in [23]; 
though subject to large uncertainties, the values are physically motivated and respect the relevant constraints (the first of which is the spectrum measured by Fermi). The 
(much lower) expected anisotropy for a continuous source distribution throughout the Galaxy (calculated with GALPROP [24]) is also shown for reference. 


single (or a few) nearby sources dominate the CRE spectrum at 
high energy. 

Key to these measurements are Fermi’s large geometric factor, 
its excellent directional capabilities and its characteristic of an 
all-sky monitor providing uniform exposure. 3 

Fermi has published the most stringent upper limits to the 
anisotropy of the CRE flux, based on the first year of data [23], 
The paper discusses ah aspects of the analysis: the pixelization of 
the sky, the pixel-to-pixel fluctuations, the search for anisotropies 
on different angular scales (from 10° to 90°) and a full spherical 
harmonic analysis. In this review we concentrate on a single as- 
pect, namely the upper limits on a possible dipole anisotropy, 
which is physically relevant for a single source dominating the 
spectrum (in which case the source defines one of the poles of 
the dipole). 

Some of the basic results are shown in Fig. 3. Though entering a 
detailed discussion of specific models and related parameter space 
is beyond the scope of this review, we emphasize that Fermi is al- 
ready starting to constrain physically motivated scenarios and 
might hope to detect a positive signal during the lifetime of the 
mission (5-10 years). If indeed there is a nearby source of high-en- 
ergy electron/positron pairs, this signal would be a powerful tool to 
shed light on its nature. For most of the dark matter scenarios one 
would expect an excess toward the galactic center (potentially at 
the same level as the ones shown in the figure) while Monogem 
and Geminga, which are two good pulsar candidates, are roughly 
in the opposite direction. 

5. Diffuse galactic gamma rays as probes of galactic cosmic rays 

Unlike the sky at visible wavelengths, the y-ray sky is strongly 
dominated by diffuse radiation originating in our Milky Way Gal- 
axy, making the galactic plane the most striking feature in any 
large-scale y-ray skymap. Although the diffuse galactic y radiation 
is brightest along the plane, it is seen at all galactic latitudes. This 
radiation is largely produced by cosmic-ray interactions with the 
interstellar gas and photon fields through the processes of inelastic 


3 The exposure is not actually uniform over the sky; there is a non-uniformity at 
the level of 1 5%, mainly due to the fact that the instrument is not taking data while 
inside the South Atlantic Anomaly. This is described in [23], 



Fig. 4. Spectral energy distribution (SED) of diffuse y radiation at intermediate 
galactic latitudes. Overlaid are SEDs for the component processes of the diffuse 
emission. The UIB component is unidentified isotropic background, including 
instrumental background. The primary component, shown by the heavier line 
above the gray band, is decay of n° mesons [25], 

nucleon scattering, bremsstrahlung, and inverse Compton scatter- 
ing. Fig. 4 shows the Spectral Energy Distribution (SED) of this 
radiation in part of the sky, along with modeled components 
[25], The largest contributor to this radiation is nucleon-nucleon 
collisions (basically cosmic ray protons and heavier nuclei hitting 
interstellar hydrogen nuclei) with subsequent decay of 7t° mesons 
into y rays in the Fermi LAT energy range. The LAT measurements 
do not confirm the unexpected excess in the few GeV energy range 
measured by EGRET on the Compton Gamma Ray Observatory com- 
pared to models. 

Mapping the galactic y rays offers a way to derive information 
about the spatial distribution of cosmic-ray particles. The modeling 
necessary to extract such results is not simple because it requires 
knowledge of all the matter and photon distributions with which 
the cosmic rays interact. Models for distribution of cosmic-ray 
sources and propagation include GALPROP [26] 4 and a derivative, 
DRAGON [27], Fig. 5 shows one result that has been derived from 
the LAT results: an indication that conventional models for CR source 


4 see also http://galprop.stanford.edu/ 
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Fig. 5. CR source distribution from one GALPROP model (solid line), compared with 
the SNR distribution obtained by the 2-D relation [31 ] (see [32] for a critique of this 
method, however) and that traced by the pulsar distribution [33] shown by dotted 
lines. The thin solid line represents an example of the modified distributions 
introduced to reproduce the emissivity gradient observed by the LAT [28], 


distributions underpredict the observations in the outer Galaxy. As 
noted in [28], the observed y-ray emissivity gradient is inconsistent 
with all available tracers of massive star formation. In this case the 
LAT results suggest either a large scale height for the CR diffusive 
halo or a flatter distribution of the sources beyond the Solar circle, 
although other possibilities cannot be ruled out, including a large 
mass of unseen gas in the outer Galaxy [29], a non-uniform diffusion 
coefficient [27], or advective CR transport by a galactic wind [30]. 

Some galaxies not exhibiting apparent nuclear jet activities are 
also seen as y-ray sources that offer information about cosmic rays 
beyond our own Galaxy. Examples are normal galaxies, including 
the Large Magellanic Cloud [35], the Small Magellanic Cloud [34], 
and starburst galaxies, such as M82 and NGC 253 [36]. GeV y rays 
in these galaxies come primarily from the interactions of cosmic 
ray hadrons and electrons with interstellar matter and photon 
fields. Fig. 6 shows the measured spectrum from the Small Magel- 
lanic Cloud. The spectrum is modeled predominantly by n° decay, 
implying hadronic processes similar to those in our own Galaxy. 
The detailed analysis of the LMC shows that the CR flux there is less 



Fig. 6. SED ofy radiation from the small magellanic cloud. Overlaid are SEDs for the 
component processes of the diffuse emission. The primary component is decay of n° 
mesons [34], 


than a third that of the Milky Way [35], and the derived CR flux in 
the SMC is even lower [34]. The y-ray luminosities of these normal 
and starburst galaxies show an approximately linear relationship 
with the product of the supernova rate and the total mass of gas 
in the galaxies [36]. The suggestion is that the y-ray production re- 
quires both target material (gas) and CR production (traced by the 
supernova rate). Although this analysis is highly simplified, the 
trend supports the expectation that the principal accelerators of 
cosmic rays are in some way related to massive-star-formation 
activity. 


6. Supernova remnants and other candidate cosmic-ray sources 

6.1. SNRs as cosmic-ray sources 

Supernova remnants (SNRs) have long been thought to be the 
main sources of galactic CRs [37], accelerating CR particles at their 
expanding shock waves through the diffusive shock acceleration 
(DSA) process [38]. In this standard paradigm, SNRs are responsible 
for the production of cosmic-ray protons and nuclei up to energies 
of ~3 x 10 15 eV, transferring ~10% of the explosion kinetic energy 
into the form of cosmic-ray energy [39]. The observations of y rays 
from SNRs have been viewed as the most promising method of 
addressing the SNR paradigm for the origin of cosmic rays, through 
the measurements of the 7t°-decay (or hadronic ) y rays [40]. De- 
spite the substantial progress of TeV y-ray observations in recent 
years [7], finding a convincing case of 7r°-decay y rays remained 
tantalizingly difficult, mainly because electromagnetic radiation 
processes involving relativistic electrons (the so-called leptonic 
components) could explain the y-ray emission as well. It has 
widely been expected that the Fermi LAT could identify the hadro- 
nic component of the y-ray emission in SNRs. Indeed the observa- 
tions with the Fermi LAT have started to provide an effective means 
to disentangle the leptonic (1C scattering and relativistic brems- 
strahlung) and hadronic (7t°-decay) y-ray components. 

Gamma-ray observations with the Fermi LAT have so far re- 
sulted in the detections of GeV y-ray emission from more than 
ten SNRs [41-48]. Fig. 7 summarizes the observed GeV luminosity 
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Fig. 7. GeV luminosity (0.1-100 GeV) in units of 10 35 erg s 1 plotted as a function of 
diameter squared for the SNRs detected with the Fermi LAT. Filled circles 
correspond to the SNRs interacting with molecular clouds. Rectangles are young 
SNRs. Pup A, middle-aged without clear molecular interactions, is shown as an open 
circle. 
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in 0.1-100 GeV as a function of the diameter squared for the 
known LAT-detected SNRs. Note that the squared diameter would 
be regarded as a reasonable indicator of SNR ages, as the age of a 
remnant in the Sedov-Taylor phase is related to its radius as 
t oc R 2 5 . The remnants interacting with molecular clouds are indi- 
cated by filled circles. They constitute the dominant class of the 
LAT-detected SNRs, and they are generally more luminous in the 
GeV band than young SNRs [49]. The observed y-ray luminosity 
of the SNRs interacting with molecular clouds is so high that one 
of the main leptonic radiation channels, namely IC scattering off 
the cosmic microwave background (CMB) and interstellar radiation 
fields, becomes unlikely. For an efficiency of ~10% for CR proton 
production and an electron-to-proton ratio of 0.01, the GeV lumi- 
nosity of the 1C emission is calculated to be L lc ~ 10 34 erg s~ ] using 
the interstellar radiation fields of the solar neighborhood. The ob- 
served luminosity well exceeds this value as shown in Fig. 7. In 
what follows, we review the results obtained with the LAT in the 
context of the SNR paradigm for CR origin. 

6.2. Young supernova remnants 

The remnants of historical supernovae are generally well stud- 
ied across the electromagnetic spectrum, making them excellent 
laboratories for studying high-energy phenomena associated with 
supernova shocks, and for testing our current understanding of the 
DSA theory. The following issues can be addressed by the Fermi LAT 
observations of young SNRs: 

(i) Acceleration efficiency — About 10% of the mechanical energy 
of supernova explosions must be transferred to the kinetic 
energy of CR protons in order to maintain the density of galactic 
CRs. Since the detailed mechanisms of how thermal particles 
are injected into the acceleration process and therefore the effi- 
ciency of particle injection remain uncertain, the current theory 
does not firmly predict the amount of CR protons in supernova 
shocks. The amount of CR electrons produced at shocks is even 
more uncertain theoretically. 


(ii) Magnetic field amplification— The DSA process at supernova 
shocks depends strongly on the generation of magnetic field 
turbulence by CR themselves ahead of the shock wave. It has 
been predicted that the turbulent magnetic field in the strong 
shocks of young SNRs can be significantly amplified by CR cur- 
rent-driven instability [50] and other instabilities. Evidence for 
the magnetic field amplification was provided by high angular 
resolution observations with the Chandra X-ray Observatory. 
Narrow widths of synchrotron X-ray filaments [51,52] and their 
short time variability [53] require strong magnetic fields. Accel- 
eration of protons up to the knee in the CR spectrum at 
~10 15 eV can be realized only if magnetic field amplification 
occurs. Gamma-ray observations can put an important con- 
straint on the strength of magnetic field, since the y-ray flux 
due to IC scattering is scaled simply by magnetic field strength 
given the observed flux of synchrotron radio emission. 

(iii) Acceleration spectrum— In the test-particle regime, the DSA 
theory predicts the spectrum of CRs accelerated at SNR shocks 
obeying a power law Q(E ) tx E~ s with s ~ 2. The CR spectra mea- 
sured directly at Earth have the form N(£) oc £~ 27 below the 
knee energy, as confirmed by the Fermi LAT measurements of 
the galactic diffuse y-ray emission. Then the diffusion coeffi- 
cient for CR propagation should be roughly D(£) oc ()(£)/ 
N(£) oc £°' 7 . Such a strong energy dependence would be incon- 
sistent with the CR anisotropy measurements [54,55]. Recently 
it has been suggested that the finite velocity of scattering cen- 
ters (usually neglected in the DSA theory) due to magnetic field 
amplification may lead to steeper spectra of accelerated CRs 
with s~ 2.2 [56], though the prediction depends on poorly 
understood properties of the magnetic turbulences produced 
at shocks. In this regard, measurements of the GeV y-ray spec- 
trum, especially that of hadronic origin, are of great importance. 

6.2.3. Cas A 

Cas A is thought to be the remnant of a Type lib supernova that 
exploded around AD 1680. Since it is at the early Sedov-Taylor 
phase, the accelerated CR spectrum can reach already its maximum 



Fig. 8. Gamma-ray spectra of young SNRs: Cas A (red) and RXJ1713. 7-3946 (black). Cas A: the LAT spectrum in 0.5-50 GeV (filled circles) [47] is shown along with the TeV 
spectra reported by HEGRA(line) [57], MAGIC (squares) [60], and VERITAS (open circles) [61]. RXJ1 713.7-3946: the LAT spectrum in 0.5-400 GeV (filled circles) [48] is shown 
in combination with the energy spectrum of H.E.S.S. -detected emission (small filled circles) [62], The y-ray spectrum of Cas A is best described by a hadronic emission model 
(7i°-decay y rays), while the y-ray emission from RXJ1713. 7-3946 can be fitted by a leptonic model (inverse-Compton scattering). (For interpretation of the references to color 
in this figure legend, the reader is referred to the web version of this article.) 
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attainable energy and the process of energy conversion, from the 
expansion kinetic energy into CR energy, can well proceed. Cas A 
was the first SNR detected in the TeV y-ray band [57]. Though 
the emission mechanisms responsible for the TeV y rays remained 
unsettled, their detections provided direct evidence for particle 
acceleration beyond 10 TeV. Observations with EGRET resulted in 
setting upper limits on the y-ray flux above 100 MeV. 

A clear detection of Cas A at a significance level of 12<r was 
made with the Fermi LAT using the y-ray data during the period be- 
tween 2008 August 4 and 2009 September 4 [47]. The GeV y-ray 
data were consistent with a point-like source, as expected from 
the angular size of the remnant (a radius of 2.5'), and also consis- 
tent with steady emission. The y-ray spectrum measured for Cas 
A in the GeV-TeV band is shown in Fig. 8. 

The X-ray observations of synchrotron-emitting filaments mea- 
sure the magnetic field at the forward shock; B ~ 0.3 mG is inferred 
from the observed width of the filaments [58], which is consistent 
with the variability timescale of the synchrotron filaments [59]. By 
combining the estimated magnetic field with the synchrotron 
spectrum in the radio bands, it was shown that the leptonic y- 
ray emission model does not well explain the LAT-detected emis- 
sion [47], Instead, the GeV-TeV y-ray emission can be ascribed 
mainly to the 7r°-decay component. The total proton content at 
the current age of the remnant amounts to W p ~ 0.4 x 10 5O erg 
for the shocked gas density of n H = 10 cm -3 , which is estimated 
from the measured dynamics of the remnant. This would be the 
first case where one can estimate the amount of CRs in SNRs with 
a reasonable accuracy. The spectral slope of accelerated CR protons 
is assumed to be s = 2.3 from the radio spectral index (Fig. 8), 
which seems compatible with the GeV-TeV spectrum. Both the 
CR content and spectrum are in reasonable agreement with the di- 
rect measurements of the CR flux and the CR anisotropy. 


6.2.2. Tycho 

Tycho’s SNR is the remnant of a Type la explosion in 1572, cur- 
rently in transition from the ejacta-dominated phase to the 
Sedov-Taylor phase. Recent detection of TeV y-rays from Tycho’s 
SNR by the VERITAS Collaboration [63] has just been followed by 
the very recent 5 a detection of GeV y-ray emission with the Fermi 
LAT [64]. The LAT spectrum is characterized by photon index 
r = 2.3 ± 0.2(stat) ± O.l(sys), which smoothly connects with the 
VERITAS spectrum. 

Given the amplified magnetic field of 0.2-0.3 mG estimated for 
shock downstream [65,58] and the constraint on the ambient 
density (n H < 0.3 cm -3 ), the y-ray spectrum of Tycho’s SNR can be 
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explained only by the 7i°-decay emission. Compared with Cas A, 
the ambient density is smaller and the radio emission is much 
weaker, so that the leptonic scenario faces a more serious diffi- 
culty. The total CR proton content at the current age amounts to 
W p ~ (0.6-1. 5) x 10 50 erg for a reasonable range of SNR parame- 
ters such as distance. The accelerated CR protons and electrons 
with number index s = 2.3 are inferred from the multiwavelength 
fitting. Such a steep spectrum may be explained by the finite veloc- 
ity of the scattering waves in the shock upstream region [66]. As 
for Cas A, the LAT results obtained for Tycho’s SNR strengthen 
the case for SNR origin of the galactic CRs. 


6.2.3. RX J1713.7-3946 

SNR RXJ1713. 7-3946 (G347.3-0.5 in the Green’s SNR list [67]) is 
well known for its bright synchrotron X-ray emission that com- 
pletely dominates the radiation output from the remnant [68,69] 
and also for its bright TeV y-ray emission that well traces the 
synchrotron X-ray map [70,62]. Since its discovery, the origin of 
the TeV y-ray emission has attracted strong interest from the 
high-energy astrophysics community [71,72]. 

Fig. 9 presents the Test Statistic (TS) map of a sky region around 
SNR RXJ1 71 3.7-3946 calculated using the events above 0.5 GeV ta- 
ken from 2008 August 4 to 2010 August 4 [48]. (TS is defined as 
TS = ~21n(L 0 /fpoint) where L 0 is the likelihood value of the null 
hypothesis that takes account of only background sources in the 
model, and Lp 0in t is that of the hypothesis assuming a point source 
at each location.) The extended y-ray emission coincident with SNR 
RXJ1713. 7-3946 can be seen in Fig. 9, though it was hard to inves- 
tigate internal spatial structures due to the lack of photons. 

In Fig. 8, the energy spectrum of SNR RX J1 71 3.7-3946 
measured with the Fermi LAT [48] is shown together with that of 
the H.E.S.S. telescopes [62]. The LAT spectrum can be characterized 
by a hard power law with photon index F = 1.5 ± O.l(stat) 
±0.1(sys), smoothly connecting with the steeper TeV spectrum. 
The hard power-law shape agrees with the expected IC spectrum 
(the leptonic model), as illustrated in Fig. 8. Given the energy flux 
ratio of the observed synchrotron X-ray emission and the y-ray 
emission, the leptonic model requires that the average magnetic 
field be weak, B~ 10 pG [62]. If the leptonic scenario is the case, 
the filamentary structures and variability in X-rays [53] would be 
attributable to locally enhanced magnetic fields, e.g., time-variable 
local turbulent magnetic fields [73]. Importantly, within the lep- 
tonic scenario the y-ray spectrum provides a robust estimate of 
the total amount of relativistic electrons, as W e ~ 1 x 10 48 erg. 
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Fig. 9. Test Statistic maps for a point source hypothesis at each location of the sky in the region around SNR RXJ1713. 7-3946 [48]. Events above 0.5 GeV are analyzed. TS 
maps were obtained by maximum likelihood fitting, where (a) the galactic and isotropic diffuse emissions and 1FGL catalog sources [74] were included as background 
sources, (b) sources A, B, and C were also taken into account in the background model. 
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Fig. 10. Fermi LAT 2-10 GeV count maps (in celestial coordinates) of the MC-interacting SNRs with extended '/-ray emission: (a) W51C; (b) W44; (c) IC443; and (d) W28 |49], 
The intensity is in units of counts per pixel with a pixel size of 0.05 x 0.05 deg 2 . Superposed are the contours from the VLA radio maps (see [49] and references therein). A 
black ellipse in panel (a) represents the location of shocked CO clumps. The black crosses in panel (c) are the locations of shocked molecular clumps from which OH maser 
emission is detected. The position of a pulsar wind nebula is also marked in panel (c). 


6.3. SNRs interacting with molecular clouds 

The brightest GeV y-ray sources associated with SNRs are mid- 
dle-aged remnants that are interacting with molecular clouds 
(MCs) as evidenced by the detections of 1720 MHz OH maser emis- 
sion [75]. Fig. 10 presents the raw LAT count maps in 2-10 GeV of 
the four SNRs (W51C [41], W44 [42], IC 443 [43], and W28 [44]) 
obtained with 2.5 yr data of the LAT observations [49]. Even before 
the subtraction of the intense galactic diffuse emission, the y-ray 
signals from the SNRs can be easily seen. 

The MC-SNRs (namely, SNRs interacting with MCs) constitute 
the dominant class of LAT-detected SNRs (Fig. 7). The y-ray lumi- 
nosity in the 1-10 GeV band spans L y = (0.8 -9) x 10 35 erg s _1 , lar- 
ger than for the young SNRs like Cas A and RX J1 713.7-3946 
(Fig. 11). The GeV y-ray spectra commonly exhibit a spectral steep- 
ening in the Fermi LAT band, demonstrating the importance of LAT 
observations of these SNRs. Also, the synchrotron radio emission of 
the four SNRs can be characterized by a large flux density of 160- 
310 Jy at 1 GHz with flat spectral index of a~0.3-0.4 [67]. The 
LAT-detected SNRs are generally radio-bright objects [49]. 

The predominance of the MC-SNRs among the SNRs detected by 
the Fermi LAT and their high y-ray luminosity indicate that the y- 
ray emission should be enhanced by the interactions with molecu- 
lar clouds. Nonthermal bremsstrahlung by relativistic electrons or 
7r°-decay y-rays produced by high-energy protons are the two 
plausible channels of the y radiation. As discussed above, the ob- 
served high luminosity of the GeV y-ray emission disfavors the IC 
scattering off the CMB and interstellar radiation fields as the emis- 
sion mechanism. 

There are two different types of models to account for the GeV 
y-ray emission from the MC-SNR systems. The Runaway CR model 
considers y-ray emission from molecular clouds illuminated by 



E [GeV] 


Fig. 11. Gamma-ray spectra of shell-type SNRs measured with the Fermi LAT: W51C 
(filled circles [41]); W44 (open circles [42]); IC 443 (filled rectangles [43]); W28 
(open rectangles [44]); Cassiopeia A (filled diamonds [47]). 


runaway CRs that have escaped from their accelerators, namely 
SNRs [76-78]. If the y-ray emission detected by the Fermi LAT 
can be ascribed to the escaping CRs, one can learn about how 
CRs are released into interstellar space and how they propagate 
in the vicinity of the SNR where self-generated Alfven waves 
may change the local diffusion coefficient [79]. Another scenario 
is the so-called Crushed Cloud model [80,81] that invokes a 
“shocked" molecular cloud into which a radiative shock (typically 
v s ~ 100 km s -1 ) is driven by the impact of the SNR's blastwave. 
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Fig. 12. Radio (left) and y-ray (right) spectra of SNR W44 together with the reacceleration model [81], 


The cosmic-ray particles accelerated at a cloud shock are adiabat- 
ically compressed behind the shock front, resulting in enhanced 
synchrotron and 7t°-decay y-ray emissions. As shown in Fig. 12, 
reacceleration of pre-existing CRs at a cloud shock is capable of 
explaining both radio synchrotron and GeV y-ray emission in the 
case of SNR W44 [81]. The y-ray luminosity of L y ~ 10 35 ergs _1 
in 1-10 GeV agrees with the theoretical expectation. Morphologies 
of synchrotron radio emission support this scenario. The radio con- 
tinuum map of SNR W44 exhibits filamentary and sheet-like struc- 
tures of synchrotron radiation well correlated with the shocked H 2 
emission [82]. One can learn about shock-acceleration processes in 
dense environments, which complements the studies of young 
SNRs. 

Both sites of the y-ray emission could be important. For exam- 
ple, a hard GeV y-ray source [44] which is also bright in the TeV 
band [83] outside the southern boundary of SNR W28 might be ex- 
plained by runaway CRs, while 7t°-decay emission from radiative- 
ly-compressed layers of clouds could account for the GeV y rays 
observed within SNR W28. Detailed analysis of spatial distribu- 
tions of GeV y rays combined with CO observations will make it 
possible to discern the y-ray sources that are attributable to run- 
away CRs. 

6.4. Other candidate cosmic-ray sources 

Although SNR remain primary candidates for accelerating the 
bulk of cosmic rays, other galactic sources observed by the Fermi 
LAT also represent sites of particle acceleration and interaction. 
Such sources include pulsar wind nebulae [84], high-mass binary 
systems such as LS 1 +61°303 [85] and Cygnus X-3 [86], whose 
spectra seem most consistent with leptonic sources, and Nova 
V407 Cygni [87], a symbiotic binary system whose spectrum can 
be modeled with either leptonic or hadronic processes. While the 
energetics requirement for the CR protons is generally difficult to 
satisfy by these non-SNR sources, they might be the dominant 
sources of the CR electrons and positrons. The large number of 
unidentified sources found in the Fermi LAT 1FGL catalog [74] of- 
fers the possibility that there may be yet-undiscovered CR produc- 
tion sources. 

7. Conclusions 

Both electron and y-ray observations made with the Fermi Large 
Area Telescope offer new insights into cosmic-ray acceleration and 


interactions. The direct electron plus positron measurements sug- 
gest, although they do not prove, a possible nearby source. Future 
LAT observations will help measure separately the positive and 
negative components up to 200 GeV and will improve on the 
search for anisotropies already begun. On a galactic scale, y rays 
suggest a larger scale height or flatter radial distribution of cosmic 
rays than previously expected. From the y-ray studies of supernova 
remnants, LAT has found growing evidence that at least some SNR 
are likely accelerators of hadronic cosmic rays, although one exam- 
ple, RX J17L3. 7-3946, appears more consistent with being a lep- 
tonic source. The Fermi LAT continues to collect data. Improved 
analysis methods, coupled with the higher statistics from the 
ongoing mission, can confidently be expected to produce even 
more results addressing questions about cosmic-ray origin and 
propagation in the future. 
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